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Abstract 

We consider the chiral expansion of the threshold amplitude for the reaction 
7P — TT^TT^p to order 0{M^). We substantiate a claim that this photoproduction 
channel is significantly enhanced close to threshold due to pion loops. A precise 
measurement of the corresponding cross sections is called for which allows to test 
chiral perturbation theory. 
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1. Single pion photo- and electroproduction off protons have been investigated extensively 
both from the experimental and the theoretical sides over the last years. Complementary 
information can be gained by the study of double pion production by real and virtual 
photons. In Ref.|]I|] we presented an exploratory study of double pion photo- and elec- 
troproduction close to threshold in the framework of heavy baryon chiral perturbation 
theory. For the threshold amplitudes, the first two terms in the chiral expansion of orders 
one and M^r, respectively, were given. The most prominent result of that investigation 
was the finding that the production cross section for final states with two neutral pions 
is considerably enhanced due to chiral (pion) loops. In this letter, we wish to sharpen 
the prediction for this effect by calculating all next-to-leading order terms of order 

in the threshold amplitude. At this order, one also has to include the effects of higher 
dimension operators. The coefficients of a subset of these have been determined recently 
from ttN scattering data in the study of the reaction vrA^ -ktiN Furthermore, one 
has to account for the Roper resonance which plays a very important role in models that 
attempt to describe the energy dependence over a large range of excitation energies 
Such studies are complementary to ours which is based on the low energy effective field 
theory of QCD and thus allows to investigate the strictures of the broken chiral symme- 
try for small external momenta, i.e. it is valid exclusively in the threshold region. The 
existing data of the two-pion photoproduction reactions from DAPHNE [§[ and TAPS 
are not yet conclusive for the threshold region. However, recently new precision data have 
been taken with TAPS for the reaction 7p — > TT^n^p also close to threshold which will 
be analysed in the near future. An improved calculation is therefore called for to allow 
new tests of chiral perturbation theory in this particular low energy reaction. 

2. First, we have to give some necessary definitions. The threshold matrix element for 
double neutral pion photoproduction off protons in the center of mass frame takes the 
form 

T(7p ^ Tr\'pf^ = la ■ (ex k) Moo , (1) 

where a, e and k denote the spin of the nucleon, the photon polarization vector and the 
photon three-momentum (in the cm frame), in order, subject to the transversality con- 
straint e-k =0. The (complex) threshold amplitude A^oo encodes the physics. Assuming 
the photoproduction amplitude to be constant in the threshold region and approximat- 
ing the three-body phase space and photon flux by a simple analytic formulap^, the 
unpolarized total cross section near threshold is 

with the photon lab energy and /i = Mt^/tti = 0.144 the ratio of the (neutral) pion 
to the proton mass. The value of the threshold photon energy in the lab frame is i?*^'' = 
2M^o(l + /i) = 308.8 MeV. 

3. To calculate the chiral expansion of the threshold amplitude A^oo^ we make use of 
heavy baryon chiral perturbation theory 0, |]. Our notation is identical to the one used 

"^■^This approximation has been shown to work excellently in 
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in Ref . [|T| and we refer to that reference for further details. We will extend the calculation 
of that paper by working out all terms of order in the Lagrangian, which amounts to 
tree graphs with insertion from the dimension one, two, three and four effective pion- 
nucleon Lagrangian as well as pion loop graphs with at most one dimension two insertion. 
Since the quantities e and k count as small momenta, we therefore get all terms up-to- 
and-including 0{M^) in the chiral expansion of A^oo- 

Let us now consider the various contributions. The tree graphs (sequential photon 
absorption and pion emissions on the proton line) with the lowest order couplings together 
with the higher order magnetic moment terms are conventionally called the Born terms. 
The Born contribution to the threshold amplitude is easily evaluated (notice that for 
double neutral pion production the first non-vanishing Born term starts at order M^r), 



M 



Born 



1-— (1 + 3.,) 



(4.9 - 1.2) GeV-^ = 3.7 GeV^^ (3) 



using eV47r = 1/137.036, = 92.4 MeV, = 139.57 MeV, m = 938.27 MeV, Kp = 
1.793 and qa = 1-32 as determined from the Goldberger-Treiman relation (with g^^N = 
13.3). Isospin breaking via the difference of the neutral and charged pion masses is 
negligible here (3%-effect). The one-pion loop diagrams contributing at order are 
taken from \M, 



= 1^1? + l + ^[2V3~ ln(2 + v^)]} = (26.9 + 16.2^ GeV-=^ 



(4) 



where the imaginary part is related to the rescattering type of diagrams which have a 
right-hand cut starting at the single pion production threshold, s = (M^ + m)^, which is 
well below the two pion production threshold. Notice that the 0{q^) loop contribution 
Eq.(^ is much larger than the one from the Born terms, Eq.(^. This is, indeed, the 
large chiral loop effect which enhances the cross section for the photoproduction of two 
neutral pions off protons considerably. At this order, one could also have 1/m suppressed 
tree graph contributions proportional to the low-energy constants q. As shown in [|l|, 
these vanish due to the threshold selection rules or cancel pairwise. However, at one order 
higher, one has such contributions. In the notation of Ref. , the terms proportional to 
the low energy constants q read 
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-2ci + (1 + Kp) 4 + (1 + 2Kp) 4' - C3 



18.4 GeV"^, (5) 



using the central values of the q as determined in p| from a fit to the subthreshold pion- 
nucleon scattering amplitudes, the ttN cr-term and some vriV P-wave scattering volumes 
(ci = —0.64, c'2 = —5.63, C2 = 7.41, C3 = —3.90, all in GeV^^). In a microscopic picture, 
these tree terms ~ q subsume all s-, t- and u-channel resonance excitations in the ir^p 
scattering amplitude. Combining these with photon absorption on the proton line, one 
is lead to a contribution to 7p —>■ tt^tt^p, i.e. the expression Eq.(^ for double neutral 
pion production off protons. The anatomy of such a chiral coefficient is discussed in M, 
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see also [T^. We note that contributions from single vr^-photoproduction via resonance 
exchanges (see Ref . |]T3[ ) followed by 7r° emission off the proton line vanish to the order we 
are working here. 

Next, we consider the pion loops at order q^. Many of these are similar in structure 
to the ones evaluated in for the reaction irN — >■ ttttN. Furthermore, one can show 
that the sum of all one-loop diagrams with one insertion proprotional to q {i = 1, 2, 3) 
vanishes. Consequently, only the dimension two operators with the LEG C4 (central value 
C4 = 2.25 GeV^^ @) together with the anomalous magnetic moments of the proton (np) 
and the neutron (k„ = —1.913) and 1/m suppressed vertices give a contribution to the q'^ 
one loop result 

-MSo"'°°' = gJ^'^4 {2[^p(ll^A - 1) + '^ni^gl + 1) + 20^7^ - 12mc4 - 4] In ^ 

~{29gl + 24mc4 + 6)— + g\7T + 3{4g\ - Amc^ - 1)^2 ln(l + V2) 

lb 

-3{gl + 2mc4 + ^) \n\l + V2) + gl{13 + llKp + 3/t„)V31n(2 + ^3) 

3 0^ 
--^^ln^(2 + V3) + 7 + 24mc4 + {g\ - l)/t„ - y (37 + 31/tp) + Hp 

(9 + llKp + 3Kn)V3 + y ln(2 + VS) I 
10.3- 17.6 i)GeV-^ (6) 
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Here, A is the scale of dimensional regularization. In what follows, we set A = m. The 
(spurious) scale dependence is rather mild, for A = m* = 1.44 GeV, the real part in Eq.(^ 
changes to —13.3 GeV~^. We note that in contrast to the calculation, the loops are no 
longer finite. Furthermore, the imaginary part at orders q^ and are of the same size 
but of opposite sign. These first two terms in the chiral expansion of Im A^oo therefore 
do not give a reliable estimate of the imaginary part. To get a better handle on it, we 
make use of the Fermi- Watson final-state theorem. Unitarity relates the imaginary part 
Im A^oo to the rescattering process 7p — > vrA^ — >• tt^tt^p with intermediate nN states on 
mass-shell. The resulting exact expression is: 



. _3M.(2 + 3,)(2 + ,) 
°° 8(l + 2/i)(l + Ai) 

3M^(2 + 3/i)(2 + /i) 



8(l + 2/i)(l + /i) 



2 

A* /Vf(l/2) /Tq a* 



7GeV-^ (7) 



where the factor in front of the brackets is the ratio of the intermediate pion momentum 
squared to the photon momentum, (in the center-of-mass system). Di^2 are the 

threshold amplitudes for the reaction tcN —>■ ttttN and the Mi_ are the pertinent single 
pion photoproduction P-wave magnetic dipole amplitudes. In the second line of Eq.(0), 
the decomposition into channels with definite total isospin 1/2 and 3/2, respectively, of 
the intermediate 71N system is given. The threshold TinN amplitudes have the (central) 
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values 



Dl = 2.26fm^ ^2 = -9.05fm^ Ao = 8.01M-^ ^32 = 2.531/"=', (8) 

and for the photoproduction magnetic multipoles we take the values from the analysis of 
Berends and Donnachie Jill at = 310.. .320 MeV (see also Ref.pl), 

pMil^'^ ~ 1 ■ 10"^ M-\ m[^1^^ ~ -6 ■ 10"=^ M-^ , (9) 

which leads to the value of Im A^oo given in Eq.(0). The imaginary parts of the (actually 
complex) T^TcN threshold and pion photoproduction magnetic dipole amplitudes are neg- 
ligible for our estimate in Eq.(7), since the respective vrA^-phases {Su ~ 2°, 5^1 — —4°) 
are tiny. Even though the imaginary part Im A^oo calculated to order is numerically 
not reliable, it serves as an important check on the evaluation of the loop integrals. 

Finally, we have to account for tree level insertions from the dimension three and four 
local contact terms. As detailed in Ref. |l|, we make use of the resonance saturation 
hypothesis to estimate these. As argued before, single A excitation (see Fig. 5a-h in 
|]I|) is to order subsumed in the tree terms ~ Cj, Eq.(|^), so we are left with double A 
graphs like in Fig. 5i,j,k in |]l|. To proceed, we generate the AA7 coupling by minimal 
substitution in the free A Lagrangian and the unknown AAvr coupling with help of the 
spin-3/2 axial current. Using the coupling constant relation (^TrAA = /5 from the 

SU(4) quark model, we get 

-^^0^ = ^^^|^v^^i(l + 2r)(l-2Z) + ^^ + f^(l-2Z)^ 

+ (1 - 2Z)ib + 2Z) + - - Z - = 1.2 GeV-=^ (10) 

3mA 4 J 

where the first term related to the graphs Fig. 5i,k in |1| (the photon does not couple 
to two A+) dominates. The off-shell parameters Y and Z have been determined in |l^ 
{gi = 5, y = 0.1, Z = —0.2). Inserting these, we find that the double A contribution is 



negligible. This agrees with the findings of Bolokhov et al. for the reaction vrA^ — > 
ttttN. We remark that such terms are exactly zero in the non-relativistic limit (static 
isobar model) 0. Another possible contribution comes from loop diagrams where the 
photon couples to a virtual meson and S° and A hyperons appear in the intermediate 
state. The respective calculation in chiral SU (3) gives after expanding in inverse powers 
of the kaon mass 

-loop 

using Mx = 493.6 MeV and D = 0.76, F = 0.50 for the axial vector coupling constants. 
As in single pion photoproduction |]l5l |l^ the (frozen) if "'"-loop effects are negligibly 



^K.-,o„p_£(C!:^F!)M = 1.4GeV- (111 



small. 

The largest resonance contribution at order comes indeed from the Roper resonance 
iV*(1440). To be specific, the pertinent Feynman diagrams involve the N^Nirn S-wave 
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vertex which is discussed in detail in Ref.0. To leading order, it is parametrized in terms 
of two coupling constants, denoted cl and Cj, respectively. We also need the N*^^^ cou- 
pling, K* , defined via the relativistic vertex (efi:*/4m)^'cr^i,F^"^\E'*. The transition magnetic 
moment k* can be determined from the Roper radiative decay width, 

V{N*+^jn) = - -1 = ^— |Ai/2|^ = O.lSMeV, K* = -0.56 (12) 



Tim" 



using the value ^41/2 = —0.072 GeV from the latest PDG tables and = (m*^ — 
m?)/2m* = 414 MeV. The strong couplings constants only known within some 

broad ranges. We therefore give an upper limit for the Roper contribution using the 
maximal value of c]; + C2 = -4.9 GeV"^ [0 

^-.af^M£i±|)<„,7GeV-. (13) 

The value Mqq = 11.7 GeV^'^ in Eq.(]TBD also includes the possible enhancement factor 

.45 which occurs if the energy denominators are not chirally 
expanded.!^ If one chooses e.g. the central values of the Roper couplings c^; + = —1.6 
GeV~^ as in this contribution is strongly reduced to J^qq = 2.6 GeV^^. Obviously, the 
Roper contribution has presently the largest theoretical uncertainty. Higher resonances 
play therefore no role for the threshold amplitude A^oo? i-e- their contribution is well 
within the theoretical uncertainty. 

4. We are now in the position to calculate the total cross section for 7p tt^tt^p 
in the threshold region. Consider first the real part of A^oo- We find that the leading 
contribution of order is at least a factor 1.5 larger than the first correction at order M^, 
Re A^oo — 31-8 GeV~'^ and for the maximal value of the Roper contribution. Re Aigo — 
21.2 GeV~'^, in order. For the central Roper couplings of pf the order correction is 
even further reduced to Re A^qo — 12.1 GeV~^, which shows that the chiral expansion is 
well-behaved, i.e. the next-to-leading order correction is smaller than the leading term. 
For the maximal iV*(1440) contribution together with the empirical imaginary part, we 
have 

A^oo ^ (53 + 7i)GeV-3 , (14) 

which shows that the imaginary part is not important for the total cross section, i.e. it 
amounts to a 2%-effect. Altogether, the predicted near threshold cross section is thenp^ 



atot(i?7) <0.91nb(^^^_^j , E'^' = 308.8MeV . (15) 

For comparison, at order q^, the prefactor in Eg . ([I5|) is 0.41 nb. We stress that the 0.91 nb 
given in Eq.(^) are the upper limit. For the central Roper couplings in 0] which lead to 

^''We are grateful to Eulogio Oset for pointing this out. 

#^We use the neutral pion mass in the phase space factor as discussed in [Q. 
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a very good description of the TrvriV threshold amphtude D2, this number is reduced to 
0.63 nb. 

In summary, we have evaluated all next-to-leading order corrections to the threshold 
amplitude for the reaction 7p vr^vr^p. The dominant contribution still comes from the 
chiral loops at leading (no n- vanishing) order, q^. The higher order corrections are in 
total not very large and increase the cross section by about a factor 1.6 to 2.2 with the 
uncertainty mainly related to the size of Roper resonance contribution. In order to test 
the prediction Eq.(|T5|) a precise measurement of these near threshold cross sections is now 
called for. 
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